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An optimized perturbation theory (OPT) at finite temper- 
ature T, which resums higher order terms in the naive per- 
turbation, is developed in O(N) </> 4 theory. It is proved that 

(i) the renormalization of the ultra-violet divergences can be 
carried out systematically in any given order of OPT, and 

(ii) the Nambu-Goldstone theorem is fulfilled for arbitrary TV 
and for any given order of OPT. The method is applied for 
the 0(4) a model to study the soft modes associated with 
the chiral transition in quantum chromodynamics. Threshold 
enhancement of the spectral functions at finite T in the scalar 
and pseudo-scalar channels is shown to be a typical signal of 
the chiral transition. 
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I. INTRODUCTION 

One of the main goals of the ultra-relativistic heavy- 
ion experiments planned at RHIC and LHC jij is to ob- 
serve the structural change of the ground state of quan- 
tum chromodynamics (QCD) at finite temperature (T), 
namely the phase transition to the quark-gluon plasma. 
The numerical simulation based on the lattice QCD is 
a powerful tool to study the static nature of this phase 
transition, in which the critical temperature and the crit- 
ical exponents are actively studied In particular, 
there exit numerical evidences that the chiral transition 
for massless two flavors is of second order, although the 
case for the real world (two light quarks + one medium- 
heavy quark) is not settled yet ||. 

If the phase transition is of second order or is close 
to it, there arises long range fluctuations in both spatial 
and temporal (real-time) directions. The latter is usu- 
ally called the soft mode and has been used as a probe 
to study phase transitions of solid states and condensed 
matters ||]. 

Despite the experimental significance of the soft modes 
in QCD, the lattice QCD simulations cannot treat such 
real-time modes in a straightforward manner. This is 
why effective theories of QCD have been adopted to study 
the time-dependent phenomena (see the reviews J5||| and 
references cited therein.) However, even in tractable ef- 
fective theories such as the linear cr-model, there exist 
subtleties at finite T. In fact, the necessity of the resum- 
mation of higher order terms in perturbative expansions 
both at high T and low T has been known for a long time 



Also, the renormalization of the ultra-violet (UV) 
divergences and related issue in resumed perturbation 
theories have been discussed in the literatures especially 
for theories with spontaneous symmetry breaking (SSB) 

a- 

Recently, we have reported our analysis of a particu- 
lar resummation method and its application to the soft 
modes in QCD |jio| . The present paper contains not only 
the detailed description of our previous analysis but also 
further investigations. 

The purpose of the present paper is twofold. Firstly, we 
will develop an improved loop- wise expansion at finite T. 
Our starting point is the optimized perturbation theory 
(OPT) (or sometimes called delta-expansion, variational 
perturbation theory etc) which is a generalization of the 
mean-field method Jll| and is known to work in various 
quantum systems |12[ . Its application to field theory at 
finite T has been considered in rcf. |l3|,[l4j for the first 
time. We will further develop the idea and prove the 
renormalizability and the Nambu-Goldstone (NG) the- 
orem in 0(N) </> 4 theory at finite T order by order in 
OPT. Our second purpose is to study the soft modes 
associated with the chiral transition in QCD by taking 
into account interactions among the soft modes (mode 
couplings). The use of OPT is essential for this purpose, 
which will be demonstrated using the 0(A) a model. 

The organization of this paper is as follows. In section 
II, we introduce a loop-wise expansion on the basis of 
OPT. The renormalization of UV divergences and the 
realization of the NG theorem in this method are also 
discussed. In section III, we will apply OPT developed 
in section II for the 0(4) a model to study the spectral 
functions of the er-meson and the 7r-meson at T ^ 0. It is 
also examined the detectability of the soft modes by the 
diphoton process a — > 2j in hot hadronic matter. Section 
IV is devoted to summary and concluding remarks. 

II. OPTIMIZED PERTURBATION AT T / 

A. Necessity of resummation at finite T 

It has been known that naive perturbations either by 
a coupling constant or by number of loops break down at 
T =/= 0, and proper resummation of higher orders is nec- 
essary M. In fact, no matter what a small dimensionless 
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coupling (say A) seems to control the perturbative expan- 
sion, the powers of T compensate the powers of A, which 
invalidates the naive expansion. 

This is easily illustrated in the 4> 4 theory: 
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Let us first consider the case /i 2 > 0. The lowest or- 
der self-energy diagram Fig.[j] (A) is 0(AT 2 ) at high T. 
However, Fig. @(B) is 0(AT 2 x ^). Further more, higher 
powers of T/fj, arise in higher loops; e.g. the n-loop di- 
agram in Figj| (C) is 0(\ n T 2n - 1 / f i 2n ~ 3 ). Thus, one 
should at least resum cactus diagrams to get sensible re- 
sults at high T 0,^5). Physics behind this resummation 
is of course the Debye screening mass in the hot plasma. 

The naive loop-expansion breaks down also for /i 2 < 0. 
The tree-level mass thq in this case is defined as 
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Here we have explicitly written the argument fj, 2 in £ for 
later use. The mass independent renormalization scheme 
with the dimensional regularization is assumed in (j^). 
Just for notational simplicity, the factor k( 4 ~") to be 
multiplied to A is omitted (n is the renormalization point 
and n is the number of dimensions). In the actual calcu- 
lations below, we take the modified minimal subtraction 
(MS) scheme. 

The c- number counter term -D/i 4 , which was not con- 
sidered in jl4| , is necessary to make the thermal effective 
potential finite. Also, it plays a crucial role for renormal- 
ization in OPT as will be shown in Sec. II. D. 

The thermal effective action T[(p 2 ] is written as the 
Euclidean functional integral |17| 
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where £(T) is the thermal expectation value of <fi. Since 
fi 2 is negative and £ 2 decreases as T increases, to 2 , be- 
comes tachyonic even below the critical temperature 
T c . Therefore, the naive loop-expansion using the tree- 
level propagator ceases to work even before the symme- 
try restoration takes place, and proper resummation of 
higher loop diagrams is necessary ||. Note that, for 
T < T c . there is no reason to believe that only the 
cactus diagrams shown in Fig.l are dominant; there ex- 
ists a three-point vertex A£<^ 3 which is not negligible for 
T~£(T). 
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where J = — dT[<p]/d(p and J^ T d 4 x = J Q 1/J dr J d 3 x. 
The "naive" loop expansion at T ^ is defined as an 
expansion by S [|18| with the tree-level mass /i 2 -I- \tp 2 /2. 

Under the naive loop-expansion with (|§|) for /i 2 > 0, 
one can completely fix the renormalization constants. 
Since ultra-violet (UV) divergences do not dependent on 
T in the naive loop-expansion JljJ, A, B, C and D are 
independent of T, and are expanded as 
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B. Resummation method 
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For theories without SSB, a systematic resummation 
method to obtain a sensible " weak-coupling" expansion 
at high T was formulated and applied to gauge theory 
and (f> 4 theory successfully |16|. 

For theories with SSB, however, loop-expansion rather 
than the weak-coupling expansion is relevant, since one 
needs to treat the thermal effective potential or the Gibbs 
free energy. We find that the optimized perturbation 
theory (OPT), which was applied to finite T system in 
p3| , pT[ , can be formulated in such a way that an improved 
loop-expansion is carried out systematically. Also, the 
method leads to a transparent renormalization procedure 
and guarantees the Nambu-Goldstone theorem order by 
order in the improved loop-expansion. 

In the following, we divide our resummation procedure 
into three steps and apply it to the 4> 4 theory. The case 
for O(N) <t> 4 theory will be discussed in Section II. E. 
Step 1: 

Start with a renormalized Lagrangian with counter terms 
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The coefficients (a;, bi,ci,di) are independent of /j, 2 , since 
we use the mass independent renormalization scheme. 
Also, the UV divergences in the symmetry broken phase 
(/i 2 < 0) can be removed by the same counter terms de- 
termined for y 2 > 

The relations of A, B, C and D with the standard 
renormalization constants are A = Z — 1, B = Z^Z — 1 
and C = Z\Z 2 — 1, where Z's are defined by <f>o = y/Z<p, 
Ao = Z\X and //q = Z^fi with suffix indicating un- 
renormalized quantities. 
Step 2: 

Rewrite the Lagrangian (JsJ) by introducing a new mass 
parameter m 2 following the idea of OPT p2[ : 



(m 2 -/, 2 ) 
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This identity should be used not only in the standard 
mass term but also in the counter terms p2| ], which is 
crucial to show the order by order renormalization in 
OPT: 
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A, B, C and D in £ c were already determined in Step 1. 

On the basis of eq.(Q), we define a "modified" loop- 
expansion in which the tree-level propagator has a mass 
m 2 + Xip 2 /2 instead of /i 2 + \ip 2 /2. Major difference be- 
tween this expansion and the naive one is the following 
assignment 



0(6°), x = 0(5). 
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The physical reason behind this assignment is the fact 
that x reflects the effect of interactions. If one adopts 
an assignment, m 2 = O(S ),x — 0(5°), the modified 
loop-expansion immediately reduces to the naive one. 

As will be shown explicitly in Sec. II. D, all the UV 
divergences in the modified loop-expansion are removed 
by the counter terms determined in the naive loop- 
expansion. 

Since (0) is simply a reorganization of the Lagrangian, 
any Green's functions (or its generating functional) calcu- 
lated in the modified loop-expansion should not depend 
on the arbitrary mass m if they are calculated in all or- 
ders. However, one needs to truncate perturbation scries 
at certain order in practice. This inevitably introduces 
explicit m dependence in Green's functions. Procedures 
to determine m are given in Step 3 below. 

To find the ground state of the system, one should look 
for the stationary point of the thermal effective potential 
V(ip 2 ) defined by 
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As mentioned above, V calculated up to L-th loops 
Vz,((/3 2 ;m) has explicit m-dependence. Thus the station- 
ary condition reads 



dV L (ip 2 ;m) 
dip 
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where derivative with respect to ip does not act on to 
by definition. Eq.([l2]) gives a stationary point of Vl for 
given to. 

One may generalize Step 2 by adding and subtracting 
ao(do<fi) 2 , ai(di4>) 2 and 70 4 [Q with o>q, ol\ and 7 being 
finite parameters to be determined by the PMS or FAC 
conditions (see Step 3). a and ot\ are especially impor- 
tant for theories with fermions at finite T and chemical 
potential [fM| . We will, however, concentrate on the sim- 
plest version («o,i = 7 = 0) in the following discussions. 



Step 3: 

The final step is to find an optimal value of to by impos- 
ing physical conditions a la Stevenson such as the 
following. 

(a) The principle of minimal sensitivity (PMS): this 
condition requires that a chosen quantity calculated 
up to L-th loops (Ol) should be stationary by the 
variation of to: 



dO L 
dm 
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(b) The criterion of the fastest apparent convergence 
(FAC): this condition requires that the perturba- 
tive corrections in Ol should be as small as possible 
for a suitable value of to. 



O l - L - n = 0, 



where n is chosen in the range, 1 < n < L. 
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The above conditions reduce to self-consistent gap 
equations whose solution determine the optimal param- 
eter to for given L. Thus to becomes a non-trivial func- 
tion of (o, X and T p6| . This together with the solu- 
tion of ( ]12| ) completely determine the thermal expecta- 
tion value £(T) = (</>) as well as the optimal parameter 
m(T). Through this self-consistent process, higher order 
terms in the naive loop expansion are resumed. 

The choice of Ol in Step 3 depends on the quantity 
one needs to improve most. To study the static nature 
of the phase transition, the thermal effective potential 
Vl(^ 2 ; to) is most relevant. Applying the PMS condition 
for Vl reads 



dV L {p> 2 \m) 
dm 
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which gives a solution to = m(ip). This can be used to 
improve the effective potential at finite T [Ol; 



(16) 



V L (<p 2 ;m)^VL(<p 2 ;m( i p)). 



Also, £(T) and m(T) are obtained by solving to- 
gether with (p"5|). In this case, the following relation 
holds: dV^imiip^/diflp^ = dV(ip 2 ; m{p)) / dip\^. 

To improve particle properties at finite T, it is more 
efficient to apply PMS or FAC conditions directly to the 
two-point functions [§7|. In ref. @, FAC with L = n = 2 
was used for the boson self-energy calculated up to two- 
loops. We will adopt a similar condition in Sec. Ill when 
we analyze spectral functions of the soft modes. 



C. UV divergence in the resumed perturbation 

We briefly mention here the reason why the renormal- 
ization in resumed perturbation is not a trivial issue. 
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In the naive perturbation theory, there arises no new 
UV divergences at T ^ because of the natural cut- 
off from the Boltzmann distribution function. Therefore, 
all the UV divergences at finite T are canceled by the 
counter terms prepared at T = 0. This statement has 
been proved in imaginary-time and real-time formalisms 
& 

On the other hand, in self-consistent methods at T ^ 0, 
the situation is not so simple since the tree-level prop- 
agators have T-dependent mass (such as m(T) in the 
above) which contains higher loop contributions through 
the self-consistent gap-equation 0. 

In fact, in most of the self-consistent methods applied 
so far (except for ref. fl4|]), the renormalization is taken 
into account "after" imposing the gap-equation. This 
procedure not only makes the renormalization non-trivial 
and hard in higher orders, but also obscures the origin 
of the UV divergences. On the contrary, in OPT ex- 
plained in the previous subsection, the renormalization is 
performed "before" imposing the gap-equation. In other 
words, the UV divergences are already removed in Step 
2, and a "finite" gap-equation is obtained from the outset 
in Step 3. 



D. Renormalization in OPT 

We now prove the order by order renormalization in 
OPT. Let us first rewrite eq. (It]) as 
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Since we use the symmetric and mass independent renor- 
malization scheme (such as the MS scheme), any Green's 
function generated by C(4>;m 2 ) can be renormalized 
solely by the coefficients A, B, C and D in C(<j>; /z 2 ). 

Suppose we make a multiple insertion of the composite 
operator (l/2)x<t> 2 to the Green's function generated by 
C(4>; m 2 ). The question is whether new divergences in- 
duced by the operator insertion are made finite only by 
the last three counter terms in (|l7|). (Note that B and 
D are already fixed in Step 1, and we do not have any 
freedom to change them.) 

The above problem is obviously related to the renor- 
malization of composite operators. In fact, the stan- 
dard method tells us that necessary counter terms 
to remove the divergences induced by the insertion of 
(l/2)x4> 2 are written as 



A 2X 2 + Ajx 



(18) 



Here Z^i is the renormalization constant for the com- 
posite operator 4> 2 , and is necessary to remove the diver- 
gence in Fig. 0(A). A2 and Ai are necessary to remove 



the overall divergences in Fig. ||(B) and in Fig. ||(C), 
respectively. 

Now, one can prove that ( |l8| ) coincides with the last 
three terms in (|l7|): 

ZZ~i - 1 = B, A 2 = D, and Ai = -2Dm 2 . (19) 

The first equation is obtained by the definition B = 
Z^Z — 1 and an identity 



Zj,2 — Z ,, 



(20) 



The overall divergence of the vacuum diagram with no 
external-legs is removed by the c-number counter term 
Dm 4 in C(<p\m 2 ). Therefore, the last two equations in 
( |l9| ) are obtained as 



Al = - { ^ 2) M = - 2Dm2 > 
2A 2 = (^) 2 [Dm*] = 2D. 



(21) 
(22) 



,^£21) is given 



For completeness, an explicit proof of 
in APPENDIX A. 

Eq.([l9|) shows clearly that all the necessary counter 
terms in OPT are supplied solely by the original La- 
grangian C{4>\[i 2 ). Let us now define {\ m2 ) 

as a renormalized n-point proper vertex with insertion 
of (l/2)x<j> 2 by j-times. (Here the external momen- 
tums are not written explicitly.) The counter terms in 
( |l8| ) together with ( [l9|) assure the finiteness of T^'^. 
Since the proper vertex can be expanded as = 
<P YaZo 1l" ' eacn coefficient 7; is also finite. This implies 
that x j ^R J) 

can be made finite order by order in OPT. 
Three comments are in order here: 

(i) Because the renormalization is already carried out 
in Step 2, one obtains finite gap-equations from 
the beginning in Step 3. Our procedure "resum- 
mation after renormalization" has several advan- 
tages over the conventional procedure "resumma- 
tion before renormalization" where UV divergences 
are hoped to be canceled after imposing the gap- 
equation. The difference between the two is promi- 
nent especially in higher order calculations. 

(ii) The decomposition (^J) should be done both in the 
mass term and the counter terms. This guarantees 
the order by order renormalization in our modi- 
fied loop-expansion. In ref. Q, the order by order 
renormalization was checked up to the two-loop or- 
der in the 4 theory at high T . Our proof shows 
that this nice feature holds in any higher orders in 
OPT. On the other hand, if one keeps the origi- 
nal counter term (1 / '2) B <\? + D/x 4 without the 
decomposition (Q), L-loop diagrams with L > M 



4 



must be taken into account to remove the UV di- 
vergences in the M-loop order (see e.g. the last 
reference in [pH|). This is an unnecessary compli- 
cation due to the inappropriate treatment of the 
counter terms. 

(iii) As far as we stay in the low energy region far below 
the Landau pole, we need not address the issue of 
the triviality of the 4> A theory p8| : Perturbative 
renormalization in OPT works in the same sense 
as that in the naive perturbation. 



E. Nambu-Goldstone theorem 

The procedure and the renormalization in OPT dis- 
cussed above do not receive modifications even if the La- 
grangian has global symmetry. For O(N) 4 theory, one 
needs to replace and (f> 2 by </> = (4>o, 4>i, ■ ■ ■, (j>N-i) and 
</> 2 respectively in all the previous formulas. In the sym- 
metry broken phase of such theory, the Nambu-Goldstone 
(NG) theorem and massless NG bosons are guaranteed 
in each order of the modified loop-expansion in OPT for 
arbitrary N. To show this, it is most convenient to start 
with the thermal effective potential V(ip 2 ). By the defini- 
tion of the effective potential, V(0 2 ) has manifest O(N) 
invariance if it is calculated in all orders. 

In OPT, V calculated up to L-th loops Vz, (^ 2 ; m) 
has also manifest O(N) invariance, because our decom- 
position (H) used in (Q) does not break O(N) invari- 
ance. Once Vl has invariance under the O(N) rotation 
(<Pi — > ipi + i9 a T^ipj), the immediate consequence is the 
standard identity: 



dV L {if 2 ;m 
difj 
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d 2 I/ L (y? 2 ;m) 
dipidtpj 



T? k <Pk, 



(23) 



with T a being the generator of the O(N) symmetry. 
Eq.(p3[) is valid for arbitrary L, m and N. 

At the stationary point where the l.h.s. of (^3|) van- 
ishes, there arises massless NG bosons for T^ k ipk ^ 0, 
since the r.h.s. of ( p3[ ) is equal to — Uy 1 ^)!^^ where 
T>ij(0) is the Matsubara propagator at zero frequency 
and momentum calculated up to L-th loops. Thus the 
existence of the NG bosons is proved independent of the 
structure of the gap-equation in Step 3. 

It is instructive here to show some unjustified approxi- 
mations which lead to the breakdown of the NG theorem. 
Many of the self-consistent methods applied so far suffer 
from these problems [g0| . 

(i) Suppose one takes into account only a part of the 
diagrams for given number of loops. Then the pion 
is no longer massless even if the symmetry is spon- 
taneously broken. Although this is a trivial point, 



sometimes such approximation is adopted in the lit- 
eratures: taking only the self-energy from the four- 
point vertex and neglecting that from the three- 
point vertex is a typical example. 

(ii) Introducing m 2 in the O(N) symmetric way as (|^) 
is a key for the NG theorem to hold in each order of 
the loop-expansion in OPT. Suppose that we make 
a general decomposition such as 



2c 2 I 2 2; \ 



(24) 



with m 2 - ^ m 2 8ij. This leads to an 0(N) non- 



invariant effective potential, and the relation (|23| ) 
is not guaranteed in any finite orders of the loop- 
expansion. For example, when the O(N) symme- 
try is spontaneously broken down to 0(N — 1), 
one may be tempted to make a decomposition 



m o (*=J=0), 



m i (*— J 7^0), and 



mfj = (i^j)- This leads to an effective potential 

Dili 1 fii m o, mi) which has only 0(N-1) 
invariance. It implies that eq.(^3|) hold only for gen- 
erators which do not mix ipo with </?«o- However, 
cq. (|2^) for those generators alone are not enough to 
prove the existence of NG bosons: In fact, (fk^o m 
the r.h.s. of ( p3] ) vanishes in the 0(N — 1) symmet- 
ric ground state, and no constraints are obtained 
for pr.i(Q). 



III. APPLICATION TO 0(4) a MODEL 

Let us apply OPT to study the soft mode associated 
with chiral transition. Our main goal is to investigate 
the spectral functions of the soft modes at finite T: 
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w,k;T) = --ImD£(«,k;T). 



Here D 1 } is the retarded correlation function 



(25) 



X>£(w,k;T) = -t / d*xe^9(t){[4>(t,x), 0(0,0)]), (26) 

where (•) denotes thermal expectation value, and 0(t, x) 
is qq(t,x.) or grygrq^t, x) in QCD. 

This spectral function was first studied using the 
Nambu-Jona-Lasinio model of QCD in the large N c limit 
The analysis shows that the scalar meson a, which 
has a large width due to the strong decay a — > 2tt, de- 
creases its mass as T increases. Eventually, a shows up 
as a sharp resonance near the critical point of chiral tran- 
sition. The detectability of such resonance was studied 
in the context of the ultra-relativistic heavy ion collisions 
p2| . Also, the spectral integrals in QCD at finite T were 
studied using the operator product expansion [p3[. 
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In the following, we will adopt a toy model (0(4) linear 
a model) to study the effect of mode couplings (interac- 
tion among the soft modes) in the one-loop level at finite 
T . This model shares common symmetry and dynam- 
ics with QCD and has been used to study the real-time 
dynamics and critical phenomena M,B5[. 



A. Determination of the parameters at T — 



The 0(4) a model reads 

t=\m? - #&\-\@? + ha 

+ l -A{d$f - - ±C($ 2 ) 2 + Dp\ (27) 

with <fi — (a, 7?) . ha is an explicit chiral symmetry break- 
ing term |3q| . A, B, C and D in the one- loop order are 

ld>TT Z £ 87H £ \K>-K A £ 

where A = — 7 + log(47r), with 7 being the Euler 
constant. 

When SSB takes place (p 2 < 0), the replacement a — > 
a + £ in eq.(]27n leads to the tree- level masses of a and 7r; 



2 , ^ fl 2 2 1 " t-2 

/i + ~f , m 07r = /x + -£ 



(29) 



The expectation value £ at T = is determined by the 
stationary condition for the standard effective potential 
dV(<p)/d(pj = 0. 

Later we will take a special FAC condition in which m 2 
deviates from /x 2 only at T ^ 0, so that the naive loop- 
expansion at T = is valid. The renormalized couplings 
/1 2 , A and h can thus be determined by the following phys- 
ical conditions in the naive loop-expansion at zero T: 

(i) On-shell condition for the pion D~ l (k 2 = to 2 ) = 
0, where — 140 MeV, and D T is the causal 
propagator for the pion in one-loop order. 

(ii) Partially conserved axial-vector current (PCAC) 
relation in one-loop; / w to 2 = h\/~Z^. Here / ff = 93 
MeV, and Z w is the finite wave function renormal- 
ization constant for the pion on its mass-shell. 

(hi) The peak position of the spectral function in the 
a channel (= mg eafc ) is taken to be 550 MeV, 750 
MeV or 1000 MeV. 

m peak _ in (iii) is consistent with recent re- 

analyses of the 7T-7T scattering phase shift p7[ . However, 
our main conclusions do not suffer qualitative change by 
other choices, 750 MeV and 1000 MeV. Instead of mP eak , 
one may take the tt-tt scattering phase shift itself as a 
condition to determine parameters pit]. However, for 



the discussions in the following, such sophistication is 
not necessary. 

We still have a freedom to choose the renormaliza- 
tion point re. Instead of trying to determine optimal 
re by the renormalization group equation for the effec- 
tive potential pSfl, we take a simple and physical condi- 
tion TOo 7r =TO 7r =140 MeV which is suitable for our later 
purpose. This choice has two advantages: (a) The one- 
loop pion self-energy S 7r (fc 2 ) vanishes at the tree- mass; 



= 0, where we have 



used the condition (i) together with = m w . (b) 

The spectral function in the a channel starts from a 
correct continuum threshold in the one-loop level. (In 
the loop-expansion, p a (Lu,0) has a physical threshold at 
lu = 2TO7r=280 MeV only if = m^.) 

Resultant parameters are summarized in TABLE I. 
The spectral functions p a and p n defined in ( p5| ) at T = 
with s = ui 2 — k 2 are shown in Fig. ||. In the 7r channel, 
there are one particle pole and a continuum starting from 
the threshold ^Jsth = toott+tooo-- \fsth is the point where 
the channel 7r + a opens. In the a channel, the spectral 
function starts from the threshold 2mo7r = 280 MeV and 
shows a broad peak centered around y/s — mP eak . The 
half width of the peak is 260 MeV, 657 MeV and 995 
MeV for m P eak = 550 MeV, 750 MeV and 1000 MeV, 
respectively. The large width of a is due to a strong a- 
2ir coupling in the linear a model. The corresponding 
cr-pole is located far from the real axis on the complex s 
plane. 



B. Application of OPT 

Now let us proceed to Step 2 in OPT and rewrite 
eq.(|27|) as 



C = \[{d$f-m^)-^( 



1 



?) 2 + \x$ 2 + ha 



— Bm 2 ^ 2 
2 v 



A 



--C($ 2 ) 2 + Dm\ 



(30) 



Since X ( = to 2 — /x 2 ) is already a one- loop order, we 
have neglected the terms proportional to Bx, D\ 2 and 
D\ which are two-loop or higher orders. 

When SSB takes place (a — > a + £), the tree-level 
masses to be used in the modified loop-expansion read 



2 2 1 " t-i 2 



A 2 
6 s 



(31) 



Since to 2 will eventually be a function of T, the tree- 
masses running in the loops are not necessary tachyonic 
at finite T contrary to the naive loop-expansion (see the 
discussion in Sec. II. A). 

The thermal effective potential V(<p; m) is calculated 
in the standard manner except for the extra terms pro- 
portional to x- The Gibbs free energy G(£; m) = V((p = 
(£, 0);to) in the one-loop level reads 
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G(fcro) 



2 C 2 



i 



4! 



64tt 2 
-T 



'Oct 



In 



'Off 



2„3/2 



3m^ In 



'Off 



2„3/2 



fire 



(32) 



d 3 k 



(2tt) s 



ln(l - e- E ° /T ) + 3 ln(l 



where 



k 2 



m; 



00- 



Although this has the similar 



structure with the standard free energy in the naive loop- 
expansion, the coefficient of the first term in the r.h.s. of 
( |32| ) is /i 2 instead of to 2 . This is because we have extra 
mass-term proportional to \ i n the one-loop level. The 
stationary point £ is obtained by 



dV(tp; to) 



dipt 



= 0. 



(33) 



Since the derivative with respect to £ does not act on 
to, this gives a solution £ as a function of T and to. By 
imposing another condition on to (Step 3), one eventually 
determines both £ and to for given T. 

At finite T, the retarded propagator has a general form 



k 2 



rn, 



0<p 



E>,k;T)' 



(34) 



with k = u> — k . The spectral function is then written 
as 



P0(w,k;T) 



1 



ImE« 



tt (fc 2 -TO 2 -Re£«) 2 + (ImI]£) 2 - (35) 



The retarded self-energy is related to the 11- 
component of the 2x2 self-energy in the real-time for- 
malism P9|]; 



ReE«(u;, k; T) = Re{E" ( Wj k) + E^ 1 (cu, k; T)} 



(36) 



ImE* ( W , k; T) = tanh( — ) Im{Sj,V> k) + E^w, k; T)}. 

Here Sj > 1 (o;, k; T) is defined as a part with explicit 
T-dependence through the Bose-Einstein distribution, 
while Ey(w,k) is the part which has only implicit T- 
dependence through m(T) and £(T). In the one- loop 
level, E^ 1 can be calculated only by the 11-component of 
the free propagator, 



it CP 

0<p\ K ' ' 1 - k 2 _ ,,,1 



mi, + te 



+ 2irn B 5{k 2 - mjL ), (37) 



with ub = [e" /T - 

One-loop diagrams in OPT for Ej, 1 are shown in Fig. ^. 
Their explicit forms are given in APPENDIX B. The NG 
theorem discussed in Sec. II. E can be explicitly checked 
by comparing eq. ( j33|) and the inverse pion-propagator at 
zero momentum [D£(0, 0; T)]^ 1 . 



C. Cancellation of T-dependent infinities 

It is instructive here to show explicitly how the UV 
divergences discussed in Sec. II. D are canceled in the one- 
loop order. The divergent part of E^(w,k;T) from dia- 
grams FigJ||(/i,i, j) reads 



A 



16tt 2 £ 
A 



(i)] 



1 



to^(T) + -to 2 ct (T) + 



16ir 2 s 

= -s?[(0 



to 2 (T) 
■ ML 



6 



(38) 



where eq.(|3l|) has been used. Namely, the terms propor- 
tional to m 2 (T) in FigJJ(ft,, i) are canceled by the counter 
term proportional to Bm 2 , while the terms proportional 
to £ 2 (T) in Fig.[|(/i, i,j) are canceled by the usual counter 
terms proportional to C. In this way, the T-dependent di- 
vergences proportional to m 2 (T) newly appeared in OPT 
is automatically canceled by the T-dependent counter 



terms obtained by the shift /i 



to 



X- The divergence 



of the free energy proportional to to 4 is also canceled by 
the last counter term in (|l7j). 

Note that the divergences proportional to \ start to 
appear from the two-loop level. They are removed by 
the counter terms proportional to \ obtained by the shift 



D. FAC condition for rn 

Since we are interested in the spectral functions in the 
one-loop level, a best way to determine to 2 is to use the 
two point function in the 7r channel. In Jl4j , a FAC condi- 
tion (|l4| ) for the two-loop self-energy at zero momentum 
(L = n = 2) was taken to obtain a gap equation for <fi 4 
theory above T c . 

The corresponding condition in our model with L = 
n = 1 reads 



E«( W = 0,O;T) 
= EiV = 0,O) 



Ei 1 (« = 0,0;T) = 0. 



(39) 



This is a condition that the one-loop correction to the 
self-energy must be as small as possible in the resumed 
perturbation theory. ( Note that ImE^(w = 0,0; T) 
vanishes identically) Unfortunately, (39) is incompati- 
ble with the condition which we adopted at T = in 
Sec. III. A to find optimal renormalization point k: 

E^(w = to 0w ,0;T = 0) = 0. (40) 

A hybrid condition which does not destroy eq.(^0|) and 
simultaneously leads to a valid gap-equation is 
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ill/ 



m 0T ,O)+Ei i (w = 0,O;T) = 0. 



(41) 



The explicit form of this FAC condition can be read from 
Eq.(B2) in APPENDIX B: 



2 2 

m = ji 



o 



(42) 



o=0 



The second (third) line is from the first (second) term in 
the l.h.s. of eg. (|4l"|) . The functions I and F are given in 
APPENDIX B (J is defined as the finite part of /.) 

At T=0, the second term in the l.h.s. of eq.(41), 
E^o^lqT = 0), vanishes by definition, and eq.(^lh for- 
mally reduces to eq.(p0[). However, we calculate (ffO|) in 
the naive loop-expansion without introducing m 2 as dis- 
cussed in Sec. III. A, while ( pl| ) is calculated with m 2 even 
at T = 0. Therefore, they are consistent only when 



m 2 (T = 0) = /j 2 



(43) 



In other words, OPT with the FAC condition (|4l| ) applied 
at T = is equivalent to the naive-loop expansion. 

In the symmetric phase at high T where £(T) ~ 0, 
eq.dH) reduces to 



m 



A 



d 3 k nB(E(m)) m 2 m 
In ■ 



(2tt) 3 E(m) 



16tt 2 



, (44) 



with E(m) = \Jm 2 + k 2 . If T 2 » to 2 , the first term in 
the r.h.s. of eq.(Q) dominates and the following solution 
is obtained: 



m 2 (T)=^ 2 + ±T 



(45) 



which implies that the Debye screening mass at high T 
can be properly taken into account. Also, both eq.(f4l|) 
and eq.(|39|) have the same solution ( j45| ) for T 2 ^> m 2 and 
are consistent with each other. For realistic values of A 
in TABLE I, the condition T 2 ^> to 2 is not well satisfied 
and one needs to solve eq. ( ffi"| ) numerically which will be 
shown in Sec. FILE. 

For intermediate values of T, eq.([4l|) can effectively 
sum not only the contributions from the diagrams in 
Fig.||(a, 6, h, i), but also from those in Fig.|](c, rf, j). 
Thus, OPT can go beyond the cactus approximation 
which sums only FigJJ(a, 6, h, i). 

Three remarks are in order here. 

(i) For sufficiently high T with fixed k, eq.(ff~4|) ceases 
to have a solution. In fact, the r.h.s. of eq.(Q) is 
always larger than the l.h.s. above a limiting tem- 
perature Ti = 500,430,420 MeV for mP a eak {T = 



0) = 500, 750, 1000 MeV, respectively. This indi- 
cates that our one-loop analysis in OPT is not suf- 
ficient for T > T\. One may try the renormaliza- 
tion group improvement by choosing e.g. k = T to 
cure this problem. However, since A in TABLE I is 
rather large, one encounters the Landau pole in the 
running coupling X(T) located at T = 440,450,490 
MeV for m? a eak (T = 0) = 500,750,1000 MeV, re- 
spectively. This again sets an upper bound of T 
beyond which the one-loop analysis in OPT is not 
reliable. 

(ii) To study the possible variations of the FAC condi- 
tion eq. (fill) , we have examined the following three 
cases, (a) Taking the high T formula eq.(|45|) in 
place of eq.([4l"|). (b) Replacing the second term 
in the l.h.s. of © by d = = TO 07r ,0;T). 
(c) Replacing the second term in the l.h.s. of ( |4l| ) 
by Ci = ReE^w = m 07r ,0;T). In the case (a), 
because of the lack of self-consistency at low T, 
m 07T {T) deviates substantially from TO 7r (T). This 
leads to an incorrect threshold for the spectral- 
function in the a channel. In the case (b), a real 
solution for to 2 is not guaranteed, because C\ is a 
complex function due to the Landau damping. In 
the case (c), a discontinuity of to 2 at certain T ap- 
pears, because C2 has a cusp structure as a function 
of uj as shown in Fig. || (C). (We will discuss this 
cusp in Sec. III. G.) Therefore, within the FAC con- 
dition for the one- loop self-energy, eq.(41) is almost 
a unique choice in the sense that it gives a smooth 
and physically acceptable solution for m 2 . 

(iii) In place of the FAC condition, one may take the 
PMS condition. However, to get a sensible gap- 
equation from the PMS condition for the thermal 
effective potential V, one needs to calculate V at 
least up to two-loops |13|. Unlike the optimized 
expansion considered in the first reference in [fl3|| , 
we have both the Hartree and the Fock diagrams 
in the two-loop order. This complicates the PMS 
analyses which will be reported elsewhere. 



E. Behavior ofm(T), m O 0(T) and £(T) 

In Fig.^(A) the tree-level masses in eq.(|3l|) and m 2 (T) 
are shown for mP* ak (T = 0) = 550 MeV. m 2 (T) is not 
tachyonic and approaches to m 2 (T) in the symmetric 
phase. This confirms that our resummation procedure 
cures the tachyon problem in Sec. II. A. 

The solid line in Fig.||(B) shows the chiral condensate 
£(T) obtained by minimizing the free energy for the case 
TO 7r (T = 0) = 140 MeV, with mP a eak (T = 0) = 550 MeV. 
£(T) decreases uniformly as T increases, which is a typ- 
ical behavior of the chiral order parameter at finite T 
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away from the chiral limit. As we approach the chiral 
limit (ft — > or equivalently — > 0), £(T) develops 
multiple solutions for given T, which could be an indica- 
tion of the first order transition. This will be discussed in 
more detail in the next subsection. The critical value of 
the quark mass m"' 1 *- below which the multiple solutions 
arise is 



/ml 



(m™V< hys -) 2 



0.08, 



(46) 



where we have used Gell-Mann-Oakes-Renner relation 
to related the pion mass with the quark mass. 



phys. 



is the physical light-quark mass corresponding 



to m*' 



140 MeV. The critical temperature for 



in), 



y TO phys. 



0.08 is T c 



170 MeV. The behavior 



of £(T) for m 7r (T = 0) = 30 MeV (just below the critical 
value m™* ) is also shown by the dashed line in Fig||(B) 
for comparison. 

Fig.|(A) shows to^ (T) for mP a eak (T = 0) = 750, 1000 
MeV. The qualitative behaviors are similar to Fig. ||(A). 
The chiral condensate £(T) for m a (T = 0) = 750, 1000 
MeV is also shown in Fig.g(B). £(T) is rather insensitive 
to the change of m£ eafe as far as m^T = 0) = 140 MeV 
is imposed. 



F. Chiral limit (h = 0) 

In the chiral limit, the FAC condition and the resultant 
gap-equation are drastically simplified and some analyti- 
cal study becomes possible. Let us carry out this analysis 
to reveal the nature of the chiral transition near the chiral 
limit. 

For h = 0, the NG theorem is satisfied for given m 2 
as shown in Sec. II. E. Therefore, as far as £ ^ (the NG 
phase), the total self-energy of the pion must vanish at 
(w,p) = (0,0): 



££(0,O;T) = 0. 



(47) 



A simultaneous solution of (^7|) and the FAC condition 



££(0, 0; T) = 0, which leads to 



m 



A 
6 



'Off 



0, and m Q(T 



A 



(48) 



The stationary condition ( |33| ) has always a solution £ = 
for ft = (the Wigner phase). The gap equation to 
determine the other solutions in the NG phase is obtained 
by substituting (|4q) into (|32 



o = /i 2 + 



A 2 



96tt 



; e In 



3k e 



d 3 k ( n B (E a ) n B {E v ) 



(2tt) 3 



E„ 



(49) 



where E a = yJk 2 + \£ 2 /3 and E« = |k|. 



The numerical solution of M9J) is given by the solid 
line in FigJ^. As can be seen from the figure, there 
are two non-vanishing solutions for £ in the range T\ = 
126 MeV <T <T 2 — 153 MeV, which is a typical behav- 
ior of the first order phase transition. For comparison, 
the case slightly away from the chiral limit is shown by 
the dashed lines in Fig.f7|. 

T\ and the behavior of £(T) for T ~ T\ can be solved 
analytically by expanding ( p9[ ) in terms of £ near £ = 0: 
Only the first term, fj, 2 , and the last two terms propor- 
tional to the Bose-Einstein distribution are relevant, and 
one obtains 



12 Ait 

£(T>T 1 )~-^=(T-T 1 ). (50) 



Ti 



The existence of the multiple solutions of the gap equa- 
tion for the 0(4) a model in the mean-field approach has 
been known for a long time |4l]] . Our analyses above con- 
firm this feature within the framework of OPT. However, 
as is discussed in the second reference of ^4j, this first 
order nature is likely to be an artifact of the mean-field 
approach, since the higher loops of massless 7r and almost 
massless a are not negligible near T\, and they could eas- 
ily change the order of the transition f42^ ]. In fact, the 
renormalization group analyses as well as the direct nu- 
merical simulation on the lattice indicate that the 0(4.) 
a model has a second order phase transition p3[ ]. 

In the following, we will go back to the the real world 
with m w (T = 0) = 140MeV, where the gap-equation has 
only one solution for given T . 



G. Spectral function at T / 

In Fig.||(A),(B), we show the spectral functions 
/w(w,0;T) for T = 50,120,145 MeV with mg eafc (T = 
0) = 550 MeV. 

In the 7r-channcl, a continuum develops for < u> < 
'"■Oct — rriQn. This originates from the induced "decay" 
by the scattering with thermal pions in the heat-bath; 
7r + 7r thcrmal _ > fj_ Because of this process, the pion ac- 
quires a width - 50MeV at T = 145 MeV, while the peak 
position does not show appreciable modification. They 
are in accordance with the Nambu-Goldstone nature of 
the pion, and are consistent with other calculations based 
on the low T expansion pi[ . 

In the er-channel for < T < 145MeV, there are two 
noticeable modification of the spectral function. One is 
the shift of the tr-peak toward the low mass region. The 
other is the sharpening of the spectral function just above 
the continuum threshold starting at u) = 2mQ^(T). 

These features are simply controlled by zeros or ap- 
proximate zeros of 



Re[D*(u,0)] 



'00 



ReS^(w,0;T), (51) 
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which appears in the denominator of the spectral func- 
tion (35). Note that the imaginary part of [D^(u, 6)] _1 
is a smooth function of u> and does not develop zeros 
above the threshold. Eq.(|5l|) is plotted in Fig.||(C). For 
T < 145 MeV, Re[D^(ui, 0)] _1 has only one zero for given 
T. This zero corresponds to an "effective" mass of the o- 
meson at finite T and roughly corresponds to the position 
of the broad peak in Fig.||(B). 

On the other hand, as T increases, the cusp in the 
low to region starts to creat an approximate zero of 
Ref-D^^O)]- 1 . (At T = 145 MeV, the cusp creates 
exact zero as shown in Fig^(C).) This is why the peak 
just above the threshold develops as T increases as shown 
in FigJ^(B). The cusp originates from the a — tt — tt cou- 
pling (the fourth diagram for the sigma self-energy in 
Fig||) and is related to the the continuum threshold by 
analyticity. The position of the cusp is exactly the point 
where the continuum starts; uj — 2mo 7l (T). 

The approximate shape of the spectral function for lo ~ 
2mo7r(T) with T ~ 145 MeV can be estimated as follows. 
The first term in the denominator of ([55|) approaches zero 
smoothly as u> — > 2mo v . 



l 0<f> 



ReZ*(u, 0) 



o. 



(52) 



On the other hand, the imaginary part of the self-energy 
is a phase space factor multiplied by a smooth and non- 
zero function /: 



Ami, 



ImS^( W ,0;T) 

Substituting (|52| ) and ( |53| ) into eq.(|35|), one finds 
p a {uj ~ 2m 07r ,0; T) 

1 

(lo - 2m 07 r)- 



(53) 



(54) 



This explains the enhancement just above the threshold 
due to the phase space factor. 

There is another explanation of the threshold enhance- 
ment. Let us start with a sum rule for the spectral func- 
tion which can be proved by the spectral decomposition 
of pa,: 



i(o;,k;T)dw 2 = {A + 1) 



(55) 



where Z = A + 1 is the wave function rcnormalization 
constant and does not depend on T. Since A=Q in the 
one-loop order, the spectral integral is unity for arbitrary 
T. This fact together with the positivity of p<j, implies 
that there is a spectral concentration near the threshold 
as mP eak (T) decreases. 



Beyond one-loop, A is divergent in perturbation the- 
ory. However, one can always define a finite and T- 
independent spectral integral as 



k; T) - p+(w, k; T = 0)] du? = 0. 



(56) 



Therefore, the same argument with the one-loop case 
holds and spectral concentration near threshold will oc- 
cur even beyond one-loop. The threshold enhancement in 
Fig.||(B), although it occurs at relatively low T, is caused 
by a combined effect of the partial restoration of chiral 
symmetry (decreasing "effective" mass) and the strong a- 
2ir coupling. In the chiral limit, the continuum threshold 
starts from u> — 0, and the enhancement occurs exactly 
at the critical temperature of chiral transition. 

Similar threshold enhancement in the it channel be- 
comes prominent just below u = rriQ a — ttiq^ for T ~ 165 
MeV. The basic mechanism of this enhancement is the 
same for the cr-case except that Re[D^(a;, O)]^ 1 is a de- 
ceasing function of T. 

The spectral functions of 7r (er) at higher temperature 
exhibits the standard behavior as expected from previ- 
ous analyses |H 33 . Shown in Fig.[| are simple a and 
tt poles and a continuum at T = 180 MeV. As T in- 
creases, these poles gradually merge into a degenerate 
(chiral symmetric) states. Due to this approximate de- 
generacy, the normal decay of a through a — > 2tt and the 
induced decay of tt through tt + tt — > a are kinematically 
forbidden at high T. This is why the width of a and tt 
vanishes. 

For sufficiently high T, the system is supposed to be in 
the deconfined phase and the decay (a, tt) — * qq must oc- 
cur. This is not taken into account in the present linear a 
model. A calculation based on the Nambu-Jona-Lasinio 
model shows, however, that there is still a chance for col- 
lective modes to survive as far as T/T c is not so far from 
unity M. 



H. Diphoton emission rate through a — > 27 

As one of the experimental candidates to see the 
threshold enhancement in the a channel, we evaluated 
the diphoton emission rate from the decay a — > 27 in 
hot hadronic matter 0| . The diphoton yield (with back 
to back kinematics) per unit space-time volume of a hot 
hadronic plasma can be written as [M 



dR. 



1 



p CT (w,k = 0;T) 



1 



(57) 



where fc M = (w,k) denotes the total four- momentum of 
the diphoton. g a is the (777 vertex at zero momentum 
and F ail {k 2 ) is a corresponding form factor. g a F ail has 
a short distant contribution from the constituent-quark 
loop and a long distant contribution from the pion-loop. 
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We took the estimates given in |47[ for these contribu- 
tions. Formula for 7r° — > 27 is obtained by a replacement 
a — > 27 in J57|). In this case, is fixed by the axial 
anomaly, and i ? 7r77 is taken from an estimate using the 
chiral quark model S] . 

The main background for the above processes is the 
thermal annihilation of pions; Tt + n~ — ► 27. Diphoton 
yield from a — ► 27 and 7r° — > 27 together with this back- 
ground at T=145 MeV are shown in Fig.^O). Threshold 
enhancement in a — » 27 process is significant only in a 
narrow region of the diphoton invariant mass and in a 
narrow region of T. Similar conclusion is drawn in other 
analysis |49(| . 

IV. SUMMARY 

In this paper, we have examined the optimized pertur- 
bation theory (OPT) in detail at finite T. For theories 
with spontaneous symmetry breaking, the loop-wise ex- 
pansion in OPT is shown to be a suitable scheme to resum 
higher order terms. 

We have shown that OPT naturally cure the two major 
problems of the naive loop-expansion, namely breakdown 
of perturbation series at high T > T c and the existence 
of tachyon poles for T <T C . 

We have also shown that OPT has several advantages 
over other resummation methods proposed so far. First 
of all, the renormalization of the UV divergences, which 
is not a trivial issue in other methods, can be carried out 
systematically in the loop-expansion in OPT. This is be- 
cause one can separate the the self-consistent procedure 
(Step 3 in Sec. II. B) from the renormalization procedure 
(Step 2 in Sec.II.B) in OPT. 

Whether the Nambu-Goldstone (NG) theorem is ful- 
filled in resummation methods has been discussed in the 
literatures. We found that the loop-expansion in OPT 
can give a clear view on this problem. The NG theorem 
is a direct consequence of the invariance of the effective 
action. Since OPT presented in this paper does not break 
the global symmetry of the effective potential in each or- 
der of the perturbation, one can prove, without much 
difficulty, that the NG theorem is fulfilled in any give 
order of the loop-expansion in OPT. 

In the latter part of this paper, we have applied OPT 
to the 0(4) a model to study the spectral functions at 
finite T. The OPT in one-loop order together with a 
FAC condition for the pion self-energy, we have success- 
fully summed not only the cactus diagrams but also other 
loop diagrams. We have demonstrated that the spectral 
function of a, which does not show a clear resonance at 
T = 0, develops a sharp enhancement near 2tt threshold 
as T approaches T c . This is due to a combined effect of 
the partial restoration of chiral symmetry and the strong 
a ~ 2ir coupling. Although it is rather difficult to ob- 
serve this enhancement in the diphoton spectrum, further 



studies will be necessary to reveal the phenomenological 
implications of this phenomena. 

The basic idea of OPT examined in detail in this paper 
will also have relevance to develop an improved perturba- 
tion theory for gauge theories in which the weak coupling 
expansion is known to break down in high orders p0| . A 
generalization of OPT, such as that discussed at Step 2 
in Sec. II. B, will be necessary for this purpose. 
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APPENDIX A: COUNTER TERMS IN OPT 



Consider the Lagrangian £(<f>;m 2 ) and define if J) (A, to 2 ) as the the renormalized n-point proper vertex with 
insertion of the composite operator cf> 2 (x)/2 by j-times. (We do not write the external momentums explicitly.) Ap- 
plying the counter terms in (^|), r^' is written by the unrenormalized proper vertex Tq 1 '^ and the renormalization 
constants (Z, Z^, A± and A 2 ) as Q 

if J ' 3 (A, to 2 ) = Z^ 2 Z- 2 j T^ j \\ ,m 2 ) + (2A 2 6 j2 + A^Sno, (Al) 

where the bare quantities with a suffix "0" are related to the renormalized quantities as Ao = Z\\, to 2 , = Z^m 2 and 
4> = y/Zcf). 

To show (ppf), consider the proper self-energy if ^ (A, to 2 ) and its derivative with respect to to 2 . After a straight- 
forward algebra using (Al), one finds 



dm 2 



-.(2,0) 



(X,m z ) = -ZZjjf'^Ao.mg) = -Z„ Z^T^'{\ m a ) 



(2,1), 



(A2) 



Since if is finite, eq.(A2) implies Z^Z^i must be finite in four dimensions. Now, in the loop-expansion with the 
A/5 scheme, Z ^ and Z^i have expansions of a form 1 + s ^aL\ a i& 1 with ai containing only the powers of 1/e. This 



fact together with the finiteness of Z^Z^i for e — > immediately leads to eq.(|20|), Z^Z^i = 1. 

To show (^l]), consider the proper vertex without external legs: if ° 3 (A,m 2 ). This quantity is a sum of all the 
one-particle irr edu cible diagrams + the c- number counter term Dm 4 . Therefore, taking a derivative with respect to 
m 2 and using (Al), one finds 



9 r (Q,Q) 
<9to 2 u 



4 JVa 



d 4 X ( 



2 (X), 



1PI + 2Dm 2 = -[if 1J (A,m 2 ) - Ai] +2Dm 2 



(A3) 



Here (-)ipi denotes the one-particle irreducible contribution. V4 denotes the Euclidean four- volume: At T ^ 0, 



Va = V3/T with V3 being the three- volume. Since if"'' is finite, eq.(A3) implies that 2Dm 2 + A± must be finite in 



four dimensions. D and Aj have expansions of the form, X^ti a z^, with a/ containing only the powers of 1/e. Thus 
one arrives at eq. (pl|) , 2Dm 2 + A\ = 0. 

To show (|22j), one needs to take derivative with respect to m 2 one more time: 



d 
dm 2 



rf t) \\,m 2 ) = z 2 zl 1 



d'xf cfy ^ ] ) 1PI + 2D = [if 2) (A,m 2 ) - 2A 2 ] + 2D. (A4) 



By the similar argument as above, D — A2 = follows from eq.(A4) 



APPENDIX B: ONE-LOOP FORMULA FOR SELF-ENERGY AT T 

Formulas corresponding to FigJj read 



k) - ffi?( W) k; T) = + F « + #) + fW] + H | fhtf) + 2F^ + F^} 

+ (- iX 2 \[4 2) + 2F a 2) + F * 3) ] + ii™ 2 - M 2 ) + counter terms 

+ (-*^) 2 [ /(3) + ^ (4) + ^ (5) ] + «0™ 2 - M 2 ) + counter terms, 



(Bl) 



(B2) 



with 



12 



r(i) 



r(2) 



j(3) = 



,2e 



,2e 



.2? 



d 4 p 



(2) 



p (3) 
F (5) 



(27r) 4 p 2 - m 2 ^ + ie (p + fc) 2 - m 2 )rj} + it ' 
<i 4 p i i 

(27r) 4 p 2 - m\ a + ie (p + k) 2 - + it ' 

rf4 p 2irn B (\p \)5(p 2 - m 2 0lf) ) 
(2^) 4 (p + fc) 2 -m 2 + ie ' 
d 4 p 



(2^) 2 n B (bo| Wbo + w|)<5(p 2 - m 2 )<5((p + k) 2 - m 2 ), 

+ (m 0cr «-> TOOtt), 



(2^) 4 

d 4 p 2irn B {\p \)8(p 2 - m^) 



(2tt) 4 (p + k) 2 - m 2 n + it 
d 4 p 



(2iryn B (\p \)n B (\p + u\)8{p 2 - m 2 a )5((p + k) 2 - mgj. 



(B3) 
(B4) 
(B5) 
(B6) 
(B7) 
(B8) 
(B9) 
(BIO) 



Here e = (4 — ra)/2, p 2 = p 2 , — p 2 , ^ 2 



The explicit forms of eq.(B3), (B4j), (B5) for k 2 > are 



k 2 and = [e"/ T - l]" 1 . 



r(i) 



r(2) 



j(3) 



16^ ( ? + 1 - l0g — 
1 P l 



16tt 2L £ 
1 ,1 



log- 



'o</> 



92 (log ^ 



1+92 1 i7F )] fOT > 4TO O0 



— 2<7 2 arctan —1 for < A; 2 < 4771 2 .,, 



167T 2 £ 



[--log- 



'Off 



'Oct 



2fc 2 



^fog^ 



(BH) 
(B12) 

(B13) 



g3 (io g ga^£;:^:o : ;[ 2 a:_g^ : ;:o : ; - 2,.)] fo r ( moa+ m ^<k 2 



2q 3 (arctan fc ~" fc %+ m °- 



arctan 



7 2 1 2 2 

fc + m o CT - m o 



2 A- 2 



„ 1 n , (2fc 2 <?3 + fc 2 -"tL+ m L)(2fc 2 93 + fc 2 +mg g -mg Ir ) 1 
^ 3 S (2fe 2 f/3-fc 2 +m^-m^)(2fc 2 1/ 3-fe 2 -m^+mg 7r )J 



for (tooct - m 07r ) 2 < k 2 < (m Qo 
for < fc 2 < (toqct — m 0w ) 2 , 



m 07r ) 2 



with 



92 



\ 



4m? 



(W> 



k 2 



- 1 



93 = 



y/\ (fc 2 + m 2 CT - m 2 ^) 2 - 4fc 2 m 2 CT | 
2fc 2 



Eq.flBJ), (|B7j), (|B8|), (gg) and flBlOl) for k = read 
_ f°° dp p 2 n B (E(m 0<t> )) 



F y ' = 



F 



(2) _ 



2vr 2 E{m 0<p ) ' 
a dp p 2 n B (E(m 0tj> )) 



1 



2^ 2 E(m 0<t> ) uj 2 -AE 2 (m 0<j> ) 



+ 6(uj 2 - Am 2 ) 



1 ^ - 4m 2 ,, oj 
-16^— ns( 2 } ' 



(■■>>) 



F 



0( W 2 - 4m 2 ) 



u 2 - 4m 2 

°^ 8^ 
dp p 2 n B (E(m 0a )) 
(2ir) 2 E(m 0<T ) 



2 /• u/ \ 

{ 



(B14) 
(B15) 
(B16) 



167ra; 



(w + £(m 0(T )) 2 - £(™ 07 r) 2 (w - S(m .)) 2 - £K) 2 

C^^^y ^( '" 2+ 1:^' ) 

for < lu 2 < (m Qa - m 07r ) 2 , (m 0a + m 0v ) 2 < u? 
for (m 0CT - m 07r ) 2 < w 2 < (m 0(T + m 07r ) 2 , 



} 



(B17) 
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F (5) 



+ (m Qa <-> m 07r ) 

' ^V(- 2 +™L--L) 2 -4r<- 2 n B ( '" 2+ i : -^' ) n 5 ( '" 2 i^' ) 

for < w 2 < (moa - m 07r ) 2 , (jn 0c r + m 07r ) 2 < 
) for (m 0c r - m 07r ) 2 < lu 2 < (m 0cr + m 07r ) 2 , 



(B18) 



where E(m) = ^/p 2 + m 2 . 
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m peak ( McV ) 


^ 2 (MeV) 


A 


h (McV a ) 


K (McV) 


T (MeV) 


550 


-284 2 


73.0 


123 3 


255 


260 


750 


-375 2 


122 


124 3 


325 


657 


1000 


-469 2 


194 


125 3 


401 


995 



TABLE I. Vacuum parameters corresponding to m£ e = 550, 750, 1000 MeV 
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(A) (B) (C) 

FIG. 1. Bubble and cactus diagrams. 






(A) (B) (C) 

FIG. 2. Diagrams which contain UV divergences as a result of the multiple insertion of (l/2)x<t> 2 - (A) corresponds to a single 
insertion with two external lines. (B) and (C) have no external lines with a single insertion and a double insertion, respectively. 
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FIG. 3. Spectral functions at T = in the n channel (A) and in the a channel (B) for m% eak = 550 MeV, 750 MeV and 1000 
MeV. 
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FIG. 4. One-loop self-energy E 11 for a and ir in the modified loop-expansion at finite T. 
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FIG. 8. Spectral function in the tt channel (A) and in the a channel (B) for T = 50, 120, 145 MeV with m p a £ak {T = 0) = 550 
MeV. The real part of (D^(cu, 0; T)) _1 as a function of ui is shown in (C). 
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FIG. 10. Diphoton yield per unit space-time volume in the back to back kinematics at T=145 MeV for mP a eak (T = 0) = 550 
MeV. 
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